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THE MATZOH BALL SOUP PROBLEM: 
A COMPLETE CHARACTERIZATION 

ROLANDO MAGNANINI AND MICHELE MARINI 


Abstract. We characterize all the solutions of the heat equation that 
have their (spatial) equipotential surfaces which do not vary with the 
time. Such solutions are either isoparametric or split in space-time. The 
result gives a final answer to a problem raised by M. S. Klamkin, ex¬ 
tended by G. Alessandrini, and that was named the Matzoh Ball Soup 
Problem by L. Zalcman. Similar results can also be drawn for a class 
of quasi-linear parabolic partial differential equations with coefficients 
which are homogeneous functions of the gradient variable. This class 
contains the (isotropic or anisotropic) evolution p-Laplace and normal¬ 
ized p-Laplace equations. 

1. Introduction 

In this paper we settle a question raised by M.S. Klamkin in m and 
extended in |A12| . 

Klamkin’s conjecture (1964). Consider the heat conduc¬ 
tion problem for a solid 12, 

ut = A u in 12 x (0,oo). 

Initially, u = 0. On the boundary u = 1. The solution to the 
problem is well-known for a sphere and, as to be expected, it 
is radially symmetric. Consequently, the equipotential sur¬ 
faces do not vary with the time (the temperature on them, of 
course, varies). It is conjectured for the boundary value prob¬ 
lem above, that the sphere is the only bounded solid having 
the property of invariant equipotential surfaces. If we allow 
unbounded solids, then another solution is the infinite right 
circular cylinder which corresponds to the spherical solution 
in two-dimensions. 

L. Zalcman |Za] included this problem in a list of questions about the 
ball and named it the Matzoh Ball Soup problem. For the case of a bounded 
solid, the conjecture was given a positive answer by G. Alessandrini m- the 
ball is the only bounded solid having the property of invariant equipotential 
surfaces. 

In |MSl| . the second author of this paper and S. Sakaguchi showed that 
to obtain the spherical symmetry of the solid in Klamkin’s setting, it is 
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enough to require that the solution has only one invariant equipotential 
surface (provided this surface is a C 1 -regular boundary of domain). In a 
subsequent series of papers, the same authors extended their result in sev¬ 
eral directions: spherical symmetry also holds for certain evolution nonlinear 
equations l |MS2i IMS4 . MS61 MS7 j); a hyperplane can be characterized as 
an invariant equipotential surface in the case of an unbounded solid that 
satisfies suitable sufficient conditions ( l.MS-li lMS5| l: for a certain Cauchy 
problem, a helicoid is a possible invariant equipotential surface I jMPrSj l: 
spheres, infinite cylinders and planes are characterized as (single) invariant 
equipotential surfaces in M 3 f jMPeSj ): similar symmetry results can also be 
proven in the sphere and the hyperbolic space ( jMS4| ). 

In [ A12 j. G. Alessandrini re-considered Klamkin’s problem for a bounded 
domain in the case in which u initially equals any function uq £ L 2 (I7) and 
is zero on <912 for all times. He discovered that either uq is a Dirichlet eigen¬ 
function or fi is a ball. A comparable result was obtained by S. Sakaguchi 
|Sa.| when a homogeneous Neumann condition is in force on <912. 

The aim of this paper is to show that Klamkin’s property of having in¬ 
variant equipotential surfaces characterizes a solution of the heat equation 
without assuming any whatsoever initial or boundary condition. This is the 
content of our main result. 

Theorem 1.1. Let 12 C M. N be a domain and let u be a solution of the heat 
equation: 

(1.1) ut = Au in LI x (0,oo). 

Assume that there exists a r > 0 such that, for every t > r, u(-,t ) is 
constant on the level surfaces of u{-,t) and Du(-,t) ^ 0 in Ll. 

Then one of the following occurrences holds: 

(i) the function ip = u(-,t ) (and hence u) is isoparametric, that is there 
exist two real-valued functions f and g such that p is a solution of the 
following system of equations: 

|H(/j | 2 = f(ip) and Atp = g(p>) in 12; 

(ii) there exist two real numbers A,// such that 

u{x , t ) = e~ Xt <j)\(x) + [i, (x, t) £ n x [r, oo), 

where 

A(f\ + A (f\ = 0 in 17; 

(Hi) there exists a real number 7 such that 

u(x, t) = 'yt + w(x), (x, t) G 17 x [r, 00 ), 

where 

Aw = 7 in 17. 

As we shall see in Remark 12.11 the assumption on the gradient of the 
solution u is only technical. Without that assumption, another (trivial) 
occurrence would be that u is constant. 

The proof of Theorem 11.11 relies on and completes those contained in |A12| 
and |Sa| : there, option (iii) and the fact that initial and boundary conditions 
are unnecessary were overlooked. 
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Isoparametric functions are well-known in the literature; accordingly, their 
level surfaces are called isoparametric surfaces and can also be characterized 
as those surfaces whose principal curvatures are all constant. The classical 
results of T. Levi-Civita [Le] and B. Segre |Se| classify isoparametric func¬ 
tions in W N by their level surfaces: they can be either concentric spheres, 
co-axial spherical cylinders (that is cartesian products of an M-dimensional 
sphere by an (N — M — l)-dimensional euclidean space, 0 < M < N — 2), 
or parallel hyperplanes (affine spaces of co-dimension 1). By using this fact, 
one can conclude fairly easily that, in Klamkin’s setting, that is when (i) of 
Theorem 11.11 holds and u is constant on Bit, then the possible shapes of a 
domain 12 can be one of the following: a ball, its exterior or a spherical annu¬ 
lus; a spherical cylinder, its exterior or a cylindrical annulus; a half-space or 
an infinite strip (see |Al2| for the simply connected case). Analogous results 
can be drawn in the case we impose on u a homogeneous Neumann condition 
(see |Sa|). 

(1.2) u v = 0 on dn x (0, oo). 

Thus, a caloric function that has invariant equipotential surfaces enjoys 
some splitting property in space-time, since it is always separable (either with 
respect to addition or to multiplication). We point out that this behaviour 
is not restricted to the case of the heat equation, since it also occurs for 
other linear evolution equations, such as the wave equation, the Schrodinger 
equation or any partial differential equation connected to the heat equation 
by some integral transform. 

More interestingly, a similar behaviour holds for the following class of 
quasi-linear evolution equations: 

(1.3) ut = Qu in 12 x (0,oo), 
where the operator 

N 

(1.4) Q u = ^2 a ij{Du) u XiXj 

i„j= 1 

is elliptic and the coefficients atj(f) are sufficiently smooth a-homogeneous 
functions of £ € R ,v , a > — 1, that is such that 

(1.5) a t ,j(a £) = cr a aij(f ) for every £ € , a > 0 and i,j = 1,... ,1V; 

this will be shown in Theorem 12.21 

Important instances of (11.311 are the evolution p-Laplace equation, 

(1.6) ut = div{\Du\ p ~ 2 Du} in 12 x (0, oo); 
the normalized evolution p-Laplace equation, 

(1.7) u t = \Du\ 2 ~ p Ayv{\Du\ p ~ 2 Du} in 12x(0,oo); 
the (anisotropic) evolution h-Laplace equation, 

(1.8) ut = AhU in 12x(0,oo); 
here, 

A hU = di v{h(Du) Dh(Du)} 
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is the so-called anisotropic /i-laplacian, where h : W N —X M is a non-negative 
convex function of class C' 2 (M Ar \ { 0 }), which is even and positively homoge¬ 
neous of degree 1 . 

In Section [3] besides a brief discussion on the possible symmetries of the 
domain 17 when initial and boundary conditions are added, we will also prove 
a classification theorem, related to equation (| 1 . 8 jh for isoparametric surfaces 
in the spirit of Levi-Civita and Segre’s result iTheorem 13.11) . 

2. The splitting property 

The proof of Theorem 11.11 and Theorem 1 2.2 1 below are similar; however, we 
prefer to give them separately to draw the reader’s attention on the different 
behaviours that depend on the parameter a in m- 

Proof of Theorem \l.l[ As originally observed in | !Allj . the assumption on the 
level surfaces of u implies that, if we set p{x) = u(x,r ), then there exist a 
number T > r and a function t/:lx [r, T) —> K, with 

(2.1) rj(s, t) = s, s E M, 
such that 

( 2 . 2 ) u(x,t) = rj(ip(x),t) 
for every (x,t) E 17 x [r, T); thus, (11.11) gives that 

r] ss (ip, t ) \Dip\ 2 + r) s (ip, t ) Aip = g t (ip, t ) in 17 x [r, T). 

By differentiating in t this identity, we obtain that ip must satisfy in 17 the 
following system of equations: 

( 2 3 \ Vs(<P,t) Aip + r) ss (ip,t) \Dip\ 2 = rj t (ip,t), 

ri s t(<p,t) Aip+ r] S3t (ip,t) \Dip\ 2 = r] tt (ip,t), 

for t € [t, T). 

Notice that the necessary smoothness of the function g can be proved 
by a standard finite difference argument: in fact, one can prove that, since 
Dp 7 ^ 0, then r] E C°°(I x [r, T)), where I = (infn <p, sup^ tp) and 7j s > 0 
on I X [t,T) (see |A111 Lemma 1], [A121 Lemma 2.1] or ] 5a] Lemma 2.1] for 
details). 

As observed in [Sa] , for the system (12.31) . it is enough to consider the 
alternative cases in which the determinant 

Z7(s, t ) = r] s r] sst - 7j st r]s S 

is zero or not zero. In fact, if D(s,t ) ^ 0 at some (s,t) E / x [r, T), then 
D ^ 0 in an open neighborhood, say U x V C / X [t, T), of (s,t): (12.3p then 
implies that 

\ D( p \ 2 = fiv) and = g{p) 

at least in a subdomain 17' of 17, and the expressions of / and g are given by 
the formulas 

, VsVtt - VstVt Vtijsst - rjttVss 

f = -, g = -; 

VsVsst VstVss VsVsst VstVss 

clearly, / and g are analytic functions. Thus u(-,t) (and hence ip) is isopara¬ 
metric in an open subdomain of 17; by using the classification result for 
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isoparametric functions by Levi-Civita and Segre and the analiticity of p, 
we can then conclude that the function ip is isoparametric in the whole fi. 
Otherwise, we have D = 0 in / x [r, T). Thus, since 

Q 2 

log (rj s ) = (rj s ) 2 (r] s rj sst - rj st r] ss ) =0 in I x [r,T), 


we have that log r] s (s,t) splits up into the sum of a function of s plus a 
function of f; ( 12 . 11 ) then implies that rj s (s,t ) only depends on t, and hence it 
is easy to conclude that 


r)(s , t) = a(t) s + b(t ), (s, t) € / x [r, T), 
for some smooth functions a and b such that 


a(r) = 1 and b{r) = 0 . 

Now, we now know that 

u{x,t ) = a(t) p(x) + 6(t) 

is a solution of <o>, thus (12.31) can be written as a linear system of equations: 

(0 a\ a'(t)p(x) -a(t)A<p = 

1 ' ' a "(t) <p( x ) ~ a '(t) ^ 

The determinant of this system must be zero, otherwise p would be constant 
(in fact, it would be that p(x) is a function of f); thus, 

a(t)a"(t) — a'(t) 2 =0 for t £ [r, oo), a(r) = 1 . 

All solutions of this problem can be written as a(t) = e _A l t_T l for A £ K 
and, by going back to (12.41) . we obtain that 

Ap(x) + A ip(x) = b'(t ) e A l f ~ r l = 7 , 

for some constant 7 . Since b(r) = 0, we have: 

1 — e _A l t_r l 

b(t) = 7 --- if A 7 ^ 0 and b(t ) =7 (t — t) if A = 0. 

A 

Therefore, we have obtained: 

u(x,t) = e _A l <_r l [<p(x) — 7 /A] + 7 /A if A 7 ^ 0, 
u{x, t) = (p(x) + 7 (f — r) if A = 0. 

In conclusion, by setting (j)\ = e Xr [p — 7 /A] and p = 7 /A for A / 0, we 
get case (ii), while setting w = p — r fT for A = 0 yields case (iii). □ 

Remark 2.1. For the sake of simplicity, in Theorem 11.11 we assumed that 
Du 7 ^ 0 in the whole f2, in order to be able to deduce the necessary regularity 
for 7 . Here, we will show how that assumption can be removed. 

The only possible obstruction to the regularity of 7 is the presence of level 
sets of p whose points are all critical. Indeed 7 is always smooth in the t 
variable and is smooth in the s variable, for every semi-regular value s of p 
( i.e. a value that is the image of at least one regular point). For every open 
subset of H of regular points for p corresponding to a semiregular value of 
p , the above theorem still holds true even if we remove the assumption on 
the gradient of the solution u. 
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We now show that, in any case, the classification holds true globally in IY 
For simplicity we assume that there exists only one critical level set. 

Let then s be a value in the range of <p such that Dpfx) = 0, for every x G 
</? - 1 (s) and set H + = {x G H : <p{x) > s}, and Q~ = {iGH : p(x) < s}. 
Clearly the above classification holds true separately in Il + and in H - . Being 
cases (i), (ii) and (iii) closed relations involving continuous functions, if the 
same case occours both in Q + and II - , then it holds in the whole domain 
IY Thus, we are left with all the cases in which in IA + and H - the solution 
assumes two different representations of those given in Theorem 11.11 

We proceed by direct ispection. If case (i) occurs in Q + or in II - , as it 
has been already remarked, then <p extends to an isoparametric function in 
an open domain containing IY 

We then have only to study the case in which instances (ii) and (iii) are 
in force in H - and Il + , respectively; by contradiction, we shall see that it is 
not possible to have a critical level set. In fact, up to sum a constant, we 
can assume without loss of generality that s / 0. As shown in ] Al2l Lemma 
2 . 2 ], the presence of a critical level set implies that Atpfx) = 0, for every 
x G 12+ FI II - . Then, according to the previous computations, we get that 
0 = Aip(x) = \ip(x), and, being A / 0, we have that ip(x ) = 0 and thence 
s = 0 , that is a contradiction. 

We now consider the case of the quasi-linear operator m ■ Notice that, 
when a = 0 , this also concerns (EH) besides the heat equation. In what 
follows, it is useful to define the generalized gradient operator: 

N 

(2.5) Qu — ^ u Xi u X j • 

i,j =1 

Theorem 2.2. Let II C M N be a domain and let u G C 1 ((0, oo); C 2 (H)) be 
a solution of equation m, 

ut = Qu in Hx (0, oo), 

where the operator Q, given in (USD, is elliptic with coefficients aij(f) that 
satisfy (11.511 . 

Assume that there exists a r > 0 such that, for every t > r, u(-,t) is 
constant on the level surfaces of u{-,t) and Du(-,t ) / 0 in II. 

Then the there exists a countable set of values ai < 02 < • • • a* < • •., such 
that 

12 = U ¥ ,_1 (K)«i+1]) 

ie N 

and, for every subdomain IT C Oi+i]) one of the following cases 

occurs: 

(i) there exist two real-valued functions f and g such that tp = u(-,t ) is a 
solution of the following system of equations: 

Qip = f(ip) and Qp = g{p) in lY; 

(ii) there exist two real numbers A,/i siLch that 

u{x , t) = [1 + A (t — r)] -1 / Q ^(.x) + p, (x, t) G lY x [r, oo), 
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with 


Q4 >A + — (j>\ = 0 in Plf 

a 

if 1; 

u(x, t ) = e ^ A ( f - r ) (j)^{x) + /i, (x, t) £Pt' x [r, oo), 

with 

Q4>\ + X(j)\ = 0 in Q', 

if a = 1 ; 

("inj t/iere exists a real number 7 such that 

u(x, t) = 7 (t — r) + u>(x), (x, t) £ Q f x [t , 00 ), 

where 

Qw = 7 in f 2 '. 

Proof. The proof runs similarly to that of Theorem 11.11 We still begin by 
setting yj(x) = u(x,r) and u(x,t ) = g(p(x),t), where 7 satisfies (12.11) . Since 
Du 7 ^ 0 equation m is uniformly parabolic, by standard parabolic regu¬ 
larity (see, for instance |LU| ), we have the necessary regularity to give sense 
to the following computations. 

By arguing as in the proof of Theorem 11.11 we obtain the system of equa¬ 
tions: 

(n Z(<P,T)Q<p + Z8(<P,T)G<p = T}t{<p,T), 

Zt{p,T)Qp + £ st (p,T)Gp = ri tt (p,T), 

where £ = (ij s ) a+1 . 

At this point, the proof is slightly different from that of Theorem 11.11 If 
there exists (s,f) such that D(s,t ) = £ £ s t — hs ft 7 ^ 0, then D(u(x,t),t) 7 ^ 0 
for x £ SI' = {x G 11 : u(x, t) £ (s — 5, s + 5)} for some 5 > 0. By setting 

~r _ f,S htt ftht — _ f,st Vt 'Htt 

1 ~ ’ 9 ~ ’ 

case (i) holds true for the function u(-,t). The same fact holds true for p, 
since u(x,t ) = rj(p(x),t) and hence p and u(x,t ) share the same level sets; 
thus, there exist functions / and g such that Gp = f(p) and Qp = g(p) in 
O'. 

Now, we define J as the closure of the complement in the image of u of 
the closure of the set {s : D(s,t ) 7 ^ 0, for some t}; J is a countable union 
of disjoint intervals. Let [a, b] C J; in fl' = 1 ([a, b]), we get 

log(0 = (0 _3 (£6rt-&&) = 0 in [a,6] x [t,T), 

£(s,r) = 1 for s £ [a,b \. 

As before, we obtain 

rj(s,t ) = aft) s + 6 (f), 

with a(r) = 1 and bfr) = 0. Proceeding as before with u(x, t ) = aft) p(x) + 
bft) gives the system: 

a'ft) p{x) — a(t) a+1 Qp = —b'ft), 
a"{t) p(x) — (a + 1) a(t) a a'ft) Qp = —b"ft). 
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The determinant of this system must be zero, otherwise p would be con¬ 
stant; thus, a must satisfy the problem 

a a+1 a" — (a + 1) a a (a') 2 = 0 in [t, T), o(t) = 1. 

The solutions of this problem are for t G [r, T) 

a(t) = [1 + A (t — r)] -1 /“ if a / 0, 
a(t) = e~ x( ' t ~ T ' ) if a = 0 , 

for some A £ M, and going back to the first equation in (12.71) gives 

Qp + — p{x) = b'(t) a(t)~ (a+1 ) = 7 if a / 0 , 
a 

Q/p + A </?(.t) = £/(£) a(f ) _1 = 7 if a = 0. 

Thus, we have that 

b(t) = {1 - [1 + A {t - r)]- 1 /*} if a/0, 

A 

1 _ e - A (*-•'■) 

bit) = 7 --- if a = 0 , 

A 

for A 7 ^ 0, 

6 (t) = 7 (t - r), 

for A = 0 and any a > 0. 

Therefore, (ii) follows when A 7 ^ 0, by setting 0 a = p—'yot/X and p = 'ya/X 
for a 7 ^ 0 and 0 a = p — 7 /A and p = 7 /a for a = 0 ; (iii) follows when A = 0 , 
by choosing w = p. □ 

3. Symmetry in space 

Theorems 11.11 and 12.21 thus give us a catalog of functions having time- 
invariant equi-potential surfaces. In particular, they inform us that the shape 
of their level surfaces can be quite arbitrary, as those of eigenfunctions are. 

This arbitrariness significantly reduces if we add special initial and bound¬ 
ary conditions. For instance, if we consider the situation in Klamikin’s con¬ 
jecture, that is if is bounded and simply connected and we require that 
the solution u of m is initially 0 on and equals a non-zero constant on 
at all tiinefQ, then the only possibility is that the level surfaces of u be 
concentric spheres, as proved in |A11] and [ 5a] . If we examine the case of a 
non-simply connected bounded domain, then the connected components of 
its boundary are spheres and must be concentric, because they are parallel. 
Thus, a spherical annulus is the only possibility. In a similar way we argue 
when is Q is unbounded, thus obtaining the remaining possible solutions: a 
(solid) spherical cylinder; a hyperspace; an infinite strip; a cylindrical annu¬ 
lus; the exterior of a ball and of a spherical cylinder. 

Actually, the presence of those initial and boundary conditions allows one 
to relax the stringent requirement that all equi-potential surfaces be time- 
invariant. In fact, to obtain the spherical symmetry of them it is enough to 
require that only one of them be time-invariant, provided it is the boundary 
of a C 1 -smooth domain compactly contained in 12 (see [MSI] ). We also 

4t can be seen that that constant can be replaced by any bounded positive function 
that bounded away from zero on dti. 
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point out that this result can be obtained under quite general regularity 
assumptions on <917, as shown in |MS7| . The case in which <917 is not bounded 
is not settled: we mention |MS3j and |MS5j for some partial results in this 
direction and |MPeS| for a characterization in 3-space. 

The possible non-regularity of <917, instead, opens up the way to non- 
spherical solutions if a homogeneous Neumann boundary condition is as¬ 
sumed on <917 for all times. As shown in |Sa| . in fact, truncated spherical 
cylinders are possible solutions, if 17 is a Lipschitz domain. At this date, no 
results are known if a homogeneous Robin boundary condition is assumed: 

u v — (tu = 0 on <917 x (0, oo). 

We notice that the situations considered in Klamkin’s conjecture and in 
[Sa j rule out the occurrence of items (ii) and (iii) in Theorem n~n Thus 
one has to deal with (i), which means that the relevant time-invariant sur¬ 
faces can only be (portions) of spheres, spherical cylinders or hyperplanes. 
Of course, the same situation reproduces for any choice of initial and/or 
boundary conditions which rule out (ii) and (iii). For instance, any time- 
independent boundary condition coupled with a non-harmonic initial condi¬ 
tion rule out (iii). 

Circumstances analogous to those of Theorem 11.11 emerge in Theoren l2.2l 
where the Laplace operator is replaced by the more general quasi-linear op¬ 
erator Q defined in m with coefficients satisfying 111.511 . In particular, 
when case (i) of Theorem 12.21 occurs, in some instances we have to deal with 
the description of anisotropic isoparametric functions and surfaces. In the 
rest of this paragraph, we shall focus on the important case of the evolution 
/i-Laplace equation (11.811 : here, it is convenient to look at the function h as 
the support function of a convex body K C belonging to the class Cl 
of bodies having a C' 2 -smooth boundary with (strictly) positive principal 
curvatures. By setting 

(3.1) H = \ K2 ' 

the system of equations that appears in item (i) of Theorem 12.21 can be 
conveniently re-written as 

(3.2) DH{Dip) ■ Dip = f(p) in 17, 

(3.3) A h ip = g{ip) in 17. 

We also notice that, since H is 2-homogeneous, by Euler’s identity (13.211 can 
be re-written as 

(3.4) 2 H(Dip) = f{ip) in 17. 

We shall say that a function ip is K-isoparametric if it is a solution of (I3.2D - 
accordingly, its level surfaces will be called K-isoparametric surfaces. 

Theorem 3.1 (AT-isoparametric functions). Let K C M. N be a convex body 
of class Cl and let h denote its support function. 

Let ip £ C 2 (f7) be a K-isoparametric function. Then its level surfaces are 
of the form. 

(3.5) Dh(§> M ) x R N - 1 ~ M , M = 0,..., N — 1. 
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Remark 3.2. Equation (13.51) should be interpreted as follows: Dh(§> M ) is an 
M-dimensional submanifold of Dh( S w_1 ) = dK and the vector-valued func¬ 
tion if : Dh( S M ) x defined by tf{Dh(v),y ) = Dh(y) + j{y), 

where j is the natural inclusion of j n M iY , defines an embedding 

and its image coincides, up to homoteties, with a level surface of p. 

The proof of the Theorem 13.11 relies on the results obtained in [GM] and 
|HLMG] , which generalize the classical ones of Levi-Civita and Segre (|Le|. 
[SeJ). In this new setting, the metric defined on the submanifolds (the level 
sets of the function <p) is given by an anisotropic (non constant) operator. 
In [GM| . the authors prove a classification theorem for hypersurfaces with 
constant anisotropic principal curvatures; the proof is mainly based on a 
Cartan-type identity which forces the A'-isoparametric surface to admit at 
most two different values for principal curvatures. 


Proof. Let p{x) be a regular value of the function p , £ the level surface 
{y E : p(y) = </?(.t)} an d let Tz(£ ) be the tangent space to £ at x. We 
also introduce the K-anisotropic Weingarten operator , 


(3.6) 


W = D 2 h 


Dp \ D 2 p 

W\) WpY 


and the K-anisotropic mean curvature , 


(3.7) 


AI h(Dp) 



Dp ■ [D 2 H(Dp)\ [D 2 p\ Dp 
\Dp\ 2 


We differentiate (13.21) and (13.41) and obtain the identities (by the square 
brackets, we denote matrices): 

[D 2 H{Dp)][D 2 p\ Dp + [D 2 p\ DH(Dp) = f\p ) Dp, 

(3.8) 2 [D 2 p\DH{Dp) = f\p)Dp, 

[D 2 H(Dp)] [D 2 p\ Dp = [D 2 p\ DH(Dp). 


After straightforward computations, from the definition (13.7p . the identi¬ 
ties m, (G3D and (13. 8p imply that 

M = g{v) - / ; (<p)/2 

vTM 


that means that M is constant on £. 

We are now going to show that M is actually the trace of the A'-anisotropic 
Weingarten operator; to do this we prove the following identity 


(3.9) 


W(x) 


D 2 H(Dp(x))D 2 p(x) 

h(Dp(x)) 


on T x (£); 


in other words, we show that the two matrices coincide as bilinear forms on 

r*(£). 
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Ill fact, (13.11) and the homogeneities of h, DH and D 2 H imply that 

[D 2 h{v)][D 2 tp\ _ {D 2 H(v)][D 2 tp\ [DH(u) ® DH(v)][D 2 <p\ _ 

\Dtp\ h(v) |V<^| h(v) 3 |V</?| 

[D 2 H(Dtp)][D 2 ip] [DH{Dtp) ® DH(Dtp)][D 2 tp\ _ 
h(Dip) h(Dtp ) 3 

_ [g^OP^ligV] 1 ,, , [flff(XV)<g>£V] 

h(Dip) 2 ^ ^ h(D<p ) 3 

where, in the last equality, we used the second identity in (1.3.81) . The desired 
formula (I3.9|i is then obtained by noticing that T X (T,) lies in the kernel of 
[DH(Dtp(x)) ® Dtp(x)\, being orthogonal to Dtp(x). 

Notice that M only depends on the geometry of the level surface; indeed, 
z/(x) = Dtp(x)/\Dtp{x)\ is the normal unit vector to £ at x and the restriction 
of — [D 2 tp(x)\/\Dtp(x)\ to T x ( S) is the shape operator of E. 

Now, we claim that there exist a relatively compact neighborhood U x C X 
of x and a number 8 > 0 such that, for any y £ U x , it holds that 

(3.10) tp(y + t DH{Dtp{y))) = tp(x + t DH(Dip(x))), 

for every 0 < r < <5. 

Without loss of generality we can assume that f(<p) = 1; indeed, since S 
is a regular level surface for tp, then f{tp) > 0; by taking = F(tp) with 
F such that (F’) 2 = /, it is easy to show that if) is another isoparametric 
function, with the same level surfaces of tp, and such that H(Dip) = 1. 

To prove our claim, we first have to show that the integral curves of 
DH(Dip) are geodesics. Let U x C S be a relatively compact neighborhood 
of x (of course, Dtp does not vanish on U x ). For every y € U x , let 7 y (r) be 
the solution of the Cauchy problem 

7 y( T ) = DH(Dtp{^ y {r))), jy(O') = y, 

and let 8 be such that, for every y € U x , remains regular on [0, d]. Then 
we have that 

<P(r/y(.T)) ~ y) = f D<p(py y (a)) • 7' (a) da = 

J 0 

[ y{?)) ■ DH(Dtp{ ly (a)))da = 2 f H{ip(p/ y (a))) da = 28, 

Jo Jo 

where we used Euler identity for H and the fact that we are assuming that 

/ = I- 

Moreover, we compute that 

ly ( T ) = [ J D 2 i7( J D^(7,(r)))][ J DV(7,(T))] 7 ;(r) = 

[D 2 H{Dtp(^ y (pr)))][D 2 tp(^ y {r))] DH^Dip^yir))) = 

= \ f(v{ly(T))) [D 2 H{Dtp( ly ( T )))} Dip(pYy(r)) = 0, 

where we used (j3.8[) and the fact that we are assuming that / = 1. 

Thus, 

7 3 / (t) = V + T DH(Dtp(y)) for 0 < r < 8, 
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and hence 


ip(y + t DH(Dip(y))) — ip(y) = 25 for 0 < r < <5, 
which means that (13.101) holds. 

Therefore, we have proved that, for every r € [0, 5), the surfaces 
T, T = {y. <p(y) = *(t)} 

are parallel with respect to the anisotropic metric induced by K: also, every 
S T has constant iv-anisotropic mean curvature. 

Thus, (13.51) follows the results |GM1 Theorem 2.1] and | OM. Theorem 1.1] 
recalled in the remark. □ 
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